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Abstract 

Population reduction is often used as a control strategy when managing infectious diseases in wildlife populations in order 
to reduce host density below a critical threshold. However, population reduction can disrupt existing social and 
demographic structures leading to changes in observed host behaviour that may result in enhanced disease transmission. 
Such effects have been observed in several disease systems, notably badgers and bovine tuberculosis. Here we characterise 
the fundamental properties of disease systems for which such effects undermine the disease control benefits of population 
reduction. By quantifying the size of response to population reduction in terms of enhanced transmission within a generic 
non-spatial model, the properties of disease systems in which such effects reduce or even reverse the disease control 
benefits of population reduction are identified. If population reduction is not sufficiently severe, then enhanced 
transmission can lead to the counter intuitive perturbation effect, whereby disease levels increase or persist where they 
would otherwise die out. Perturbation effects are largest for systems with low levels of disease, e.g. low levels of endemicity 
or emerging disease. Analysis of a stochastic spatial meta-population model of demography and disease dynamics leads 
to qualitatively similar conclusions. Moreover, enhanced transmission itself is found to arise as an emergent property of 
density dependent dispersal in such systems. This spatial analysis also shows that, below some threshold, population 
reduction can rapidly increase the area affected by disease, potentially expanding risks to sympatric species. Our results 
suggest that the impact of population reduction on social and demographic structures is likely to undermine disease 
control in many systems, and in severe cases leads to the perturbation effect. Social and demographic mechanisms that 
enhance transmission following population reduction should therefore be routinely considered when designing control 
programmes. 
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Introduction 

The relevance of ecology to understanding the dynamics and 
persistence of infectious disease has long been recognised [1], and 
ecological factors are critical to wildlife disease systems. Control of 
disease in wildlife is of considerable importance for managing risks 
to humans [2,3] and livestock [4,5], as well as for the conservation 
of wildlife species themselves [3,6-8]. Population reduction is a 
commonly employed strategy used to control disease in wildlife 
[9, 1 0] with the aim of reducing the number of infected animals 
and the overall size of key populations, leading to a reduction in 
rates of transmission, disease prevalence and risks to other 
populations. Application of this strategy is supported by theoretical 
evidence of a threshold for disease persistence below which disease 
does not spread quickly enough to persist, and eventually dies out 
[9, 1 1 , 1 2] . However, there is growing evidence that population 
reduction may be less effective than standard analyses predict, and 
in some cases be counter-productive (see below). Such unexpected 
increases in disease prevalence following population reduction 
have been termed the "perturbation effect" [13]. The theoretical 
basis and empirical evidence for disease thresholds in wildlife has 



been reviewed [14], concluding that important elements of wildlife 
ecology are neglected by current theories. 

It is known that the social and spatial structure of host 
populations has significant implications for disease persistence and 
prevalence [15,16]. Population reduction disrupts existing social 
structures and this may lead to increased numbers of contacts [17] 
and/or a greater proportion of agonistic encounters within or 
between groups [18,19]. Similarly, a change in susceptibility of 
individual hosts may also occur as a consequence of population 
reduction due to stress [20]. Both effects will enhance disease 
transmission and are likely to be widespread and reduce or even 
reverse the efficacy of population reduction measures. 

For example, management of rabies in foxes ( Vulpes vulpes) has 
shown that vaccination is more suitable than culling, as the latter 
can destabilise social structure and lead to enhanced transmission 
rates [10,21]. Studies of the management of classical swine fever 
(CSF) in wild boar (Sus scrofa) recommend that hunting should 
cease following detection of the disease [22], in order to discourage 
dispersal of infected individuals, and reduce risks to neighbouring 
groups [10]. The U.K. Randomised Badger Control Trial (RBCT) 
[23] showed that reactive culling of badgers (Meles meles) in 
response to a confirmed bovine tuberculosis (Mycobacterium bovis, 
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bTB) herd breakdown in catde, was associated with a 27% 
increase in the incidence of confirmed breakdowns, relative to 
survey-only trials [24]. Repeated reactive culling was also 
associated with increased bTB prevalence in badgers [25]. 

In this paper we study the potential for behavioural and 
demographic aspects of the ecology of wildlife species to reduce or 
reverse the efficacy of population reduction as a means of disease 
control. Our results are based on the analytical and numerical 
treatment of generic models of demography and disease dynamics 
in wildlife populations. In a non-spatial context we analyse the 
potential that individual and collective behavioural responses to 
population reduction have on disease control. We use this 
framework to explore the demographic and epidemiological 
characteristics of wildlife disease systems that make them 
susceptible to such effects. We then demonstrate that such impacts 
arise as an emergent property of spatial models of wildlife disease 
systems with density dependent dispersal. Finally we discuss the 
significance of these results for disease control in wildlife. 

Methods 

A non-spatial deterministic model of demography and 
disease dynamics 

We examine a generic single pathogen wildlife disease system 
with a fluctuating host population. The number of susceptible and 
infected individuals in the population at time t are S(i) and I(t) 
respectively, and the total population size is given by Mt) — S(t)+I(f). 
We assume density dependent (logistic) growth, with intrinsic 
reproduction rate r (the maximum rate that individuals can 
reproduce in optimal circumstances), limited by a carrying 
capacity c (the population size for which the density limited per- 
capita birth rate reaches zero — note this is not necessarily the 
same as the population equilibrium [26,27], since mortality, 
including that induced by disease and population reduction, will 
prevent the population from attaining this maximum). Natural 
mortality (from causes unrelated to disease or explicit population 
reduction measures) occurs at constant per-capita rate d, while 
disease induced mortality occurs at constant per-capita rate e. The 
rate of infection is a combination of susceptibility and contact rates 
between susceptible and infective individuals and here we consider 
density dependent infection (i.e. disease transmission depends on 
the density of infectives, I) with horizontal transmission rate /?. 

We model population reduction as a constant per-capita death 
rate p which applies to all individuals regardless of disease status. 
As noted earlier such measures can alter host behaviour and hence 
contact rates. We therefore model the horizontal disease 
transmission rate as fi+kp. Here k>0 represents any mechanism 
or combination of mechanisms that lead to increased contact rates 
or susceptibility in a host population subjected to population 
reduction at rate p. Note that this formulation represents a 
simplification in that the effect is linear in p, there is no lag as p 
changes and the effect is constant for the duration of the 
population reduction event. 

In Appendix SI (see Supporting Information, File SI), we show 
how to formulate a simple non-spatial deterministic model that 
encapsulates the above assumptions. We also simplify this 
representation, removing the variables c and r by respectively 
scaling the variables S, I and M by 1 1 c to obtain values between 0 
(empty) and 1 (at carrying capacity) and rescaling time by r (see 
Appendix SI. 2, Eqn. SI in File SI). Analysis can then focus on the 
effects of population characteristics (parameters d and e), disease 
dynamics (/?), population reduction (p), enhanced transmission (k) 
and the interactions between them. However, results for specific 
values of c and r can still be obtained by appropriate back scaling. 



The rescaled deterministic ordinary differential equations (ODEs) 
that combine the demography and disease dynamics described 
above with population reduction and a corresponding enhanced 
transmission resulting from explicit behavioural and implicit 
ecological (system) responses are given by: 

S = N(l-N)-(d+p)S-(P + kp)SI (1) 



i=-(d + e+p)I + (p + kp)SI 

Three fixed points of this system of equations are derived in the 
Supporting Information: population extinction, where {S,I} = 
{0,0}; the disease free equilibrium, where {S,I} = {1 — d— p,Q}; 
and the endemic equilibrium {S,I} = {S*(p),I*(p)}, where both 
the population and the disease persist (note that it is possible for 
I*(p) to be negative, in which case 1(f)— >0, since there cannot be a 
negative number of individuals). The stability properties of these 
equilibria are discussed in Appendix SI. 4 in File SI. Note that we 
write the endemic equilibrium as a function of the reduction rate p, 
even though it also depends on other parameters, because we are 
particularly interested in the effect of population reduction. 

A spatial stochastic model of demography and disease 
dynamics 

In the stochastic spatial model we consider a set of sites where, 
at time /, the integer number of susceptibles and infectives in site i 
are S,(t) and /,(?) respectively. Since we are dealing with numbers 
of individuals these are not rescaled as above. The demography 
and disease dynamics of each sub-population are governed by the 
same processes as for the non-spatial model, with the addition of 
dispersal and disease transmission within and between groups. 

Dispersal is the movement of individuals between social groups, 
for the purposes of obtaining more resources such as food or 
reproductive opportunities (including inbreeding avoidance). In 
the model dispersal from any given site occurs at constant per- 
capita rate m, into any of its nearest neighbouring sites. However, 
since this process may be mediated by the population levels in the 
destination site [28-30] this is modified by a function f(Nj), where 
Nj is the population at neighbouring site j. We consider a step 
function 



m)= 



1 if Nj<aN^ ¥ 
0 if Nj>aN^ F 



(2) 



where = c(\ — d / r) is the population size in the disease free 
equilibrium, and a. is the fraction of the disease free equilibrium at 
which the neighbouring site becomes accessible. Dispersal rates 
may also be affected by conditions in the source area, e.g. due to 
overpopulation, social exclusion, or lack of resources, lack of 
mating opportunities in small populations; however, we do not 
consider these effects here. 

Disease transmission rates within and between groups are 
denoted fi w and /?4 respectively. The horizontal disease transmis- 
sion rate in site i is therefore given by 



where the sum is over neighbouring sites of i. The total infection 
rate is given by H= ^2, Hi and the effective disease transmission 
rate is defined as 



PLOS ONE | www.plosone.org 



2 



May 2014 | Volume 9 | Issue 5 | e86563 



Demographic Processes Counter Disease Control 



,elf (E $)(£/,) 

i i 

The spatial model is implemented as a discrete state-space Markov 
process, to account for demographic stochasticity, with events and 
associated rates shown in Table 1, and simulated using the 
Gillespie algorithm [31]. In the spatial model population reduction 
is parametrised by the probability that a site is targeted pi and the 
rate of removal of individuals within targeted sites p2- 

Measuring the perturbation effect 

We define the magnitude of the perturbation effect at time / 
after the application of population reduction at rate p to be 

n{t;p) = I(t;p)-I(t;Q) (3) 

A population that is in equilibrium /*(0) prior to the application 
of population reduction at rate p, will reach a new equilibrium 
I*(p). We define the persistent perturbation effect 

n eqm = max{T (p)fi) - max{T (0),0} (4) 



(disease endemic before intervention starts), is subjected to 
sustained population reduction at rate p (see Fig. 1A). The long 
term equilibrium I*(p) increases with k (i.e. the effectiveness of 
population reduction reduces) and when k is greater than some 
critical value hp, n eqm >0. However, another behaviour is also 
apparent: when k approaches a lower threshold k h there is a 
temporary increase in I(t), which results in Yl(t;p)>0 for a short 
period, despite no perturbation effect in the long term (Il eqm < 0). 
We call these two increases the persistent and the transient 
perturbation effect, and examine their properties in the following 
sections. Both persistent and transient perturbation effects are also 
possible in the case of emerging disease (when starting from close 
to the disease free equilibrium) (see Fig. IB). 

Behaviour in the long-term equilibrium can be seen by plotting 
the endemic equilibrium I*(p) versus population reduction rate p, 
for several values of the horizontal transmission rate p. Three 
important points are evident (see Fig. 1C). First, persistent 
population reduction at a sufficiendy intense rate does reduce 
the level of disease, leading to I*(p)<I*(0). Second, the maximum 
size of the persistent perturbation effect reduces as the horizontal 
transmission rate increases, with no perturbation effect present in 
the deterministic model for ft sufficiendy high. Finally, increased 
horizontal transmission induced by population reduction can allow 
the disease to persist, where it would otherwise fade out in the 
absence of culling. 



Note that I*(p) may be negative, in which case the equilibrium 
is no longer stable and cannot be reached, and so I(t)—*0, hence 
the restrictions (see Appendix S1.4 in File SI for more details). In 
the results, we study both the persistent II eqm and transient n(f ; p) 
perturbation effects. In the spatial case we also examine the 
proportion of sites containing infectives, Pi(t; p), as the basis for 
measuring the perturbation effect 

n shes (t;p) = P I (t;p)-P I (t;0) (5) 



Persistent perturbation effect with no disease induced 
mortality 

We now explore the properties of the persistent perturbation 
effect n eqm in more detail. For clarity we focus on the algebraically 
simpler case where there is no disease induced mortality, e = 0, and 
technical details of the analysis are given in Appendix SI. 5 (in File 
SI). Subsequently we apply numerical analysis to Eqn. 1 with 
disease induced mortality e>0. 

Case 1: Disease persists without population reduction, I (0)>0 

In this case there is a perturbation effect if 



Results 

Explicit enhancement of disease transmission induced by 
population reduction 

We first consider the perturbation effect in the deterministic 
non-spatial model. Several features of the perturbation effect 
caused by increased horizontal disease transmission in response to 
population reduction are demonstrated in Fig. 1. For different 
levels of transmission enhancement k, a range of outcomes are 
possible when a population in the endemic equilibrium 7*(0) 

Table 1. Default event rates for the stochastic SI model. 



„ „ d+p d „ 

n >= n ^=- p -WFp + e >0 

Minimum disease enhancement required to produce 
perturbation effect. Note that when k=0, we obtain 
FIi = —p(l + 1/yS), which is always negative, showing that culling 
reduces disease when there is no mechanism enhancing disease 
transmission. Rearranging gives a threshold value of k, above 
which a perturbation effect is possible 





Event 


Rate 


<fSi 


SI, 


SSj 


SI, 


Birth of S, 


rN t (\-N,/c)5t 


+1 


0 


0 


0 


Death of S, 


dSt&t 


-1 


0 


0 


0 


Death of /, 


(d + e)I,6t 


0 


-1 


0 


0 


Infection of S, 


H,St 


-1 


+1 


0 


0 


Dispersal of S, to site j 


mzSif(Nj)M 


-1 


0 


+1 


0 


Dispersal of /, to site j 


mzIJ(Nj)&t 


0 


-1 


0 


+1 



Event rates and corresponding effects in the spatial stochastic model. H, and f(A/y) are defined in the methods. 
doi:1 0.1 371 /journal.pone.0086563.t001 
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Table 2. List of parameters used in the deterministic and stochastic SI models. 





Parameter 


Symbol 


Non-spatial 


Spatial 


Intrinsic reproduction rate 


r 


1 


1 


Carrying capacity 


c 


1 


20 


Natural mortality rate 


d 


0.2 


0.01 


Disease induced mortality rate 


e 


0.1 


0.1 


Horizontal transmission rate 


P 


0.4 




background 


ft 




0 


within group 


ftv 




0.5 


between groups 


ft. 




0 


Dispersal rate 


m 




0.1 


threshold value 


a 




0.7 


Population reduction rate 


P 


0.1 




coverage 


Pi 




0.2 


removal within sites 


Pi 




0.5 


Disease enhancement 


k 


5 





A summary of the parameters and their symbols used in the non-spatial and spatial models are described here. Values shown indicate both the parameters and their 
default values used in the spatial and non-spatial models. 
doi:1 0.1 371 /journal.pone.0086563.t002 



k>k] 



1(1 +/P 

d-pji 



There is a lower bound on this threshold, such that 
k\> fi(\+fi)/d>{\+ f})> \ (since the disease is able to persist, 
which requires that P>d/(l—d)>d, hence P/d>\). 

High disease prevalence precludes a perturbation 
effect. As /?— >co, III— »— p, showing that for sufficiently high 
P, the perturbation effect cannot occur in the deterministic model. 
In fact in this case there is an upper bound, /?„, on the value of ft 
for which 111 > 0, 



1 



P<P U = ^ {\+pk? +4dk- l -{\+pk) 

and ITi >0 only when /?</?„ (see Fig. 2 A for high ff). Similarly, as 
d—*Q, Yl\ -* —p( 1 + 1 /(fi + kp)) < 0, showing that the perturbation 
effect is possible only for higher mortality rates. There is a 
corresponding lower bound on d for which ITi > 0, at 



d>d, = p(l+P+kp)/k 

(see Fig. 2B, for low d). For low d<di, infectives are removed from 
the population slowly, and for high /?>/?„, the disease spreads 
quickly; either situation leads to disease saturation, with insuffi- 
cient susceptibles to allow for a perturbation effect. 

High rates of population reduction will reduce disease 
levels. A simple observation is that a persistent perturbation 
effect is possible (for any model) only if the population size under 
persistent culling is greater than the equilibrium number of 
infected individuals without population reduction i.e. 
N*(p)> I*(0) which implies that there is an upper bound on the 
culling rate, p = d/P, above which population reduction will 
reduce disease (see Appendix SI. 5 in File SI). This is also evident 
in k\ (the lower bound for k) which diverges as />—►£?/ p from below 
implying that, in order to see a perturbation effect, population 
reduction must produce ever greater enhanced transmission k as p 



approaches this critical level. Furthermore (see Appendix SI. 5 in 
File SI), we show that the range of/) that permits the perturbation 
effect also depends on k and is given by 



o< P < Tj 



1 + 0 



This is illustrated in Fig. 1C, where the range of p for which 
I*(0)<I*(p) decreases with and that for sufficiently large p, 
I*(p) <0 for all fi. 

Case 2: Disease does not persist without population reduction, 
7*(0)>0 

Population reduction can allow disease to persist where it 
would naturally fade out. In this case there is a persistent 
perturbation effect if the disease is only able to persist under 
continued population reduction for a given p and k, i.e. when 



n 2 = n e , 



\-d- 



l + kp 



>0 



The conditions for which the disease is only able to persist under 
population reduction are detailed in Appendix SI. 4 in File SI. 
The minimum k in order to make the endemic equilibrium 
I*(p)> 0 following population reduction is 



k>k 2 



d+p-p(l-d-p) 
pQ-d-p) 



and therefore sufficientiy large k can lead to a perturbation effect 
under these conditions. For example, given d = 0.2,/) = 0, 7*(0) < 0 
for /?<0.25; however, given = 0.2,/? = 0.2, (thus unable to persist 
for p = 0), when population reduction is applied at rate p = 0.l, 
then I*(p = 0.1)>0 for k>ki =2.286, therefore the disease can 
persist as long as population reduction is sustained, leading to a 
perturbation effect (see Fig. 1C, for /? = 0.2). 
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Time (years) 



persistent increase 



control 



transient increase 



no increase 
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k=3 



k=2 
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Time (years) 



80 
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0.5 



0.4 



0.3 



I °- 2 



« 0.1 




-0.1 



-0.2 





0.2 




0.3 




0.5 




0.8 


J3u 





0.1 0.2 0.3 0.4 

Population reduction rate, p 



Figure 1. Deterministic simulation of /(f), and algebraic solution of /*(p). The results of ongoing population reduction are shown for various 
levels of disease enhancement k in (A) endemic disease, (B) emergent disease (starting near the disease free equilibrium, 1(0) = 0.05). (C) shows the 
endemic equilibrium for varying ft. The lines cut the vertical axis at /*(0), and so the perturbation effect occurs whenever a line rises above this value. 
Note that for ft = 0.2, the equilibrium is negative for small p (which cannot be reached, since only a non-negative number of individuals is biologically 
possible), and so if any disease is introduced for p = 0, it moves to the disease free equilibrium /p F = 0, and the perturbation effect does not occur 
until p is sufficiently high. The dotted line shows fi u (see text for details), marking the upper bound of ft for given p for which the perturbation effect is 
possible, and crosses each line at the point where the increase no longer occurs for that value of /J (f} u is also illustrated in Fig. 2A). Parameters are 
given in Table 2, except p = 0.2 in (A) and (B). 
doi:1 0.1 371 /journal.pone.0086563.g001 



Persistent perturbation effect with disease induced 
mortality 

We now investigate the persistent perturbation effect in the 
more complex situation with disease induced mortality e>0 by 
solving Eqn. 1 numerically to show how n eqm varies with e itself, 
and also with horizontal disease transmission background 
mortality d, and enhanced transmission k resulting from popula- 
tion reduction. 

Numerical analysis of the role of transmission rate /J is consistent 
with the analysis of the previous section (see Fig. 2 A). Under case 1 
(where /*(0)>0 and n = IIi), II decreases with /? and no 
perturbation is possible for /J > /?„ because the disease has 



saturated the population, whereas in case 2 (where 7*(0)<0 and 
II = n.2), II increases with /J, and there is a lower limit below 
which the disease becomes extinct despite enhanced disease 
transmission. This is in accordance with analysis of Yl\ and IT2 (see 
Appendix SI. 5 in File SI, and above). 

The role of natural mortality d is also consistent with the 
previous analysis (see Fig. 2B). In the region of case 1, there is a 
lower bound d/, below which the perturbation effect is not possible 
due to disease saturation, and above which Yl increases with d. In 
the region of case 2, II decreases with d, and there is an upper limit 
on d, above which the disease becomes unable to persist despite 
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0 0.1 0.2 0.3 0.4 0.5 0.6 1 10 100 1000 



Disease induced mortality, e Disease enhancement, k 

Figure 2. Sensitivity analysis of the persistent and transient perturbation effects in the deterministic model. Parameter values as in 
Table 2 and marked by a red dot when explicitly varied, the transient perturbation (green) is shown for f = 5, the persistent perturbation (blue) n cqm is 
evaluated at t= 1,000. The transient perturbation: (A) has an optimum for intermediate ji, decays with (B) natural and (C) disease induced mortality, 
and (D) increases with k. With the exception of (D), the behaviour of the persistent perturbation is more complex. In (A) to the left of ["(p) = 0 (dashed 
vertical line) ft is low and there is no perturbation (I*(p)< 0). To the right of /*(/?) = /*(0) (dotted vertical line corresponding to upper bound ji u , see 
text) the prevalence in the absence of culling is sufficiently high to prevent a perturbation. The central region between 7*(/;) = 0 and /*(/)) = /*(0) is 
divided by a third vertical line /*(0) = 0 (dot dashed), independent of k and p, into regions corresponding to case 1 (/*(0)>0) where the disease 
persists the absence of population reduction, and case 2 (7*(0)<0) where it does not (see text for details). The maximum persistent perturbation 
occurs at this boundary. Under case 2, population reduction is sufficient to stabilise the endemic equilibrium. In (B) as natural mortality d increases 
from zero (moving left to right) 7*(0) decreases, and the pattern seen in (A) is reversed. Here the dotted vertical line l*(p) = /*(0) denotes the lower 
bound d\. (C) shows the impact of disease induced mortality e is similar to that of natural mortality, but the chosen parameter values mean that 
prevalence is never too high to prevent a perturbation effect. Note: dotted and dashed lines are reversed when k is too low for the perturbation effect 
to occur, leaving no room for cases 1 and 2. See Fig. 4 for analogous spatial model results. 
doi:1 0.1 371 /journal.pone.0086563.g002 



enhanced transmission. The role of disease induced mortality e, is 
broadly similar to that of d (see Fig. 2C). 

The impact of the disease enhancement parameter k on the 
perturbation effect is illustrated for case 1 in Fig. 2D. FI increases 
with k, tending to an asymptote as fc— >oo, while there is no 
perturbation effect below the threshold k\ . The behaviour under 
case 2 (not shown) is broadly similar with a different lower bound 
kj and lower asymptote. 



Maximising the persistent perturbation effect 

An important addendum to these results is related to the 
conditions that maximise the perturbation effect. For low mortality 
rates d or e, or high transmission rate /?, the disease is able to persist 
before and during population reduction, the prevalence is very 
high and there is little room for further increase. As mortality 
increases or transmission decreases, the size of the perturbation 
effect n eqm increases until P = (d + e)/(\ — d), where the endemic 
equilibrium becomes negative, and the disease becomes unable to 
persist for p = 0 (as in case 2). After this point, as mortality 
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increases, or transmission decreases, n e q m decreases, until 
mortality is too high, or transmission is too low to maintain the 
disease either before or during population reduction. This implies 
that the maximum perturbation effect occurs when /* (0) x 0 and 
l*(p)>0. Therefore in practice, the persistent perturbation effect 
is most likely in a disease with very low prevalence. These results 
can be seen graphically in Fig. 2. 

Transient perturbation effect 

The transient perturbation effect can be assessed by linearising 
the system and examining the rate of change of Yl(t;p) with 
respect to time, at time t=0, which is positive (i.e. the disease 
increases faster under population reduction) only if Ie(0,N — l/k) 
(see Appendix S2.1 in File SI). To obtain an initial increase in 
disease levels there must be some infectives, but similar to the 
results in the persistent case, too many infectives will prevent a 
transient perturbation effect; as k increases a transient perturbation 
effect is possible for ever larger numbers of infectives. The lower 
bound here is equivalent to k>k t = \ / S and since Se(0, 1), this 
requires that k,>l; therefore, the transient perturbation effect 
does not occur in the absence of a change in behaviour. It is also 
possible to show that the transient perturbation effect increases 
fastest when S = N/2+\/2k and I = N/2-\/2k (i.e. roughly 
equal numbers of susceptibles and infectives) and that 11(0; p) 
increases with both p and k (see Appendix S2 in File S 1 for details). 
Also, a temporary peak, where I(i)> I*(p) may occur, if the 
disease increases quickly before culling reduces the population size 
N; this can be observed in both endemic and emerging disease 
cases (see Fig. 1). 

Starting from the endemic equilibrium. Consider the case 
where the disease is in the endemic equilibrium {5*,/*} prior to 
disease intervention (as shown in Fig. 1A). We show in Appendix 
S2 (in File SI) that k t = ji/(d + e) so that the minimum disease 
enhancement required for a transient perturbation effect is 
reduced when the infection rate /? is small and mortality rates d 
and e are large. In addition k, < k p (where k p is the relevant k\ or 
ki) and the transient perturbation effect occurs for smaller k than 
the persistent perturbation effect. Consequently, for small k<k t , 
there is no perturbation effect. For larger ke(k t ,k p ), I(t; p) > I(t; 0) 
for small t, i.e. the number of infectives is initially larger following 
disease intervention, however eventually I(t)-*I*(p) which is less 
than the initial level 7*(0), and in this case the increase is 
temporary. However, for k>k p , the number of infectives increases 
and remains higher than the control. 

Starting from near the disease free equilibrium. The 
situation is somewhat different in the case of an emerging outbreak 
where 1(0) = e where e > 0 is small, and S(0) = <Sq F — e, as shown 
in Fig. IB (Sq F is S in the disease free equilibrium, see Appendix 
SI. 3 in File SI). Here, k, = l/(l-d-e) (see Appendix S2.2 in File 
SI), and so the minimum disease enhancement required for a 
transient perturbation effect is reduced when the initial prevalence 
is low (although contrary to the persistent perturbation effect, 
when mortality rates are also low). Fig. 2 shows the impact of 
varying d, e, fS and k on the transient perturbation effect for the 
case of an emerging outbreak where 1(0) = e where e > 0 is small, 
and 5(0) = Spp — e. These numerical results show that the 
transient perturbation monotonically decreases with both natural 
and disease induced mortality, whilst it monotonically increases 
with enhanced transmission k. The disease transmission rate /? 
affects the time disease takes to reach equilibrium, and therefore 
small /? can result in a slow initial increase (and small transient 
perturbation effect), while very large /? can saturate the population 
and prevent the transient perturbation effect from occurring at the 
time considered; the largest increase therefore occurs with an 



intermediate value of /?, although this will vary depending on the 
time at which the transient perturbation effect is assessed. 

These results contrast with those for the persistent perturbation 
effect (also shown in Fig. 2), demonstrating that conditions 
required for the transient and persistent perturbation effect are 
not necessarily the same for both emerging and endemic disease. 

Implicit enhancement of disease transmission induced by 
population reduction 

We now show how the intrinsic dynamics of a natural spatial 
formulation of disease transmission and demography may give rise 
to an increased effective horizontal transmission when population 
reduction is applied, leading to an implicit perturbation effect. The 
non-spatial results of the previous section suggest that perturbation 
is strongest when disease prevalence is relatively low and where 
population reduction is intermediate, and gives rise to a sufficiendy 
large increase in the horizontal transmission rate. 

We begin by demonstrating the importance of heterogeneity in 
the model, and show through analysis of the horizontal infection 
rate in a simple two-site model that in the spatial model such an 
enhancement of the transmission rate will be strongest in situations 
where infection levels are most heterogeneous between groups. 

Heterogeneity and the perturbation effect in the spatial 
model 

Consider a simple two-site model, with density dependent 
dispersal between the two groups A and B. The global infection 
rate H is 

H = H A +H B = fi w (S A I A + S b Ib) + MSaIb + S B I A ) 

Assuming disease induced mortality rate e = 0, then N A = Nb = N, 
and the number of infectives is I a +Ib = I, thus H can be 
simplified to 

H = constant + 2(fi w - p h )(I - I A )I A 

revealing that when between-group infection rate pb is small, and 
Pw > fib> H i s maximised when the infection is distributed evenly 
between sub-populations, and I A =Ib- Conversely, //is minimised 
when I a =0 or I A =1 (i.e. all infectives are restricted to one of the 
groups). 

In order to quantify enhanced transmission resulting from 
population reduction, consider the rate of change of the global 
infection rate H (differentiated with respect to time) to obtain 

H = p w (S A I A + sj A + S B I B + s B i B ) 
+ MSaIb + sJ b + s b i a + s B i A ) 

While H is affected by all processes, including birth, death, 
infection and dispersal, if we only examine the effect of dispersal 
on H by substituting only the relevant components 
(S A = ... +mS B f'(N)-mS A f(N) etc.) then we obtain 

//d Isp cr S al = 2m(P w - p h )(I A - I B ) 2 f(N) 

Consequendy, if P w > P h and I A ¥=Ib, then the effect of dispersal is 
to increase H. Moreover, this rate of increase in horizontal disease 
transmission is greatest for heterogeneously distributed disease 
(larger difference I a —Ib), and for larger dispersal rate m. The 
presence of the density dependence function ffN) shows that it is 
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Figure 3. Time trajectories and heterogeneity for emergent disease in the stochastic model. (A) Population numbers, S(f), /(f), and N(f). (8) 
Proportion of sub-populations containing infectives, P/(t), effective transmission rate /J, and dispersal rate. (C) Distribution of / across sites. (D) 
Effective transmission rate /i for disease transmission vs population reduction coverage p,. Parameters are given in Table 2, and initial conditions are 
at the disease free equilibrium {5,/} = {20,0}, while in 20% of sites randomly chosen, a single individual is infected, resulting in {S,7} = {19,1}. 
Population reduction occurs annually from years 50-69, and in p\ =20% of sites (chosen randomly each year) the removal rate is set to pi = \», 
without regard to disease status (equivalent to an overall culling rate of p = 0.2). An initial reduction in / is rapidly replaced by an increase, which is 
due to the increased chance of invasion of naive groups by infectives due to the density dependent dispersal. The CI for the effective transmission 
rate increases for large p\ due to the increasing number of simulations where the disease becomes extinct. 
doi:1 0.1 371 /journal.pone.0086563.g003 



greater for smaller JV, which will follow as a consequence of 
population reduction. Note that the formula for H shows that a 
large Pi, will lead to rapid spread between sites, quickly spreading 
to sub-populations, and reducing spatial heterogeneity in the 
distribution of disease. 

Initial conditions 

Given the above discussion, when studying the spatial model we 
focus on cases where disease is distributed heterogeneously 
between groups and overall prevalence is low. This is most easily 
achieved when the system is close to the disease free equilibrium 
with: (i) disease maintained in each site by high within-site 
transmission rate ji m and low mortality; (ii) low levels of disease 
transmission between sites; and (iii) relatively large and stable 



populations at each site leading to low levels of dispersal between 
sites. Under this scenario, even in the absence of population 
reduction, the number of sites infected, and thus overall 
prevalence, tends to slowly increase (from close to the disease free 
equilibrium) as rare dispersal or transmission events spread 
disease. Fig. 3 (discussed in detail below) shows how transient 
perturbation effects occur in such a system. In contrast, we show in 
Appendix S3. 2 (in File SI) that by making both disease and 
population less stable within sites it is possible to achieve a 
dynamic quasi-equilibrium (quasi- because the ultimate fate of all 
simulations of this model is total extinction) where the spread of 
disease to uninfected sites is balanced by spontaneous recovery of 
infected sites, e.g. through death of infectives and birth of 
susceptible individuals. When the system is in such an endemic 
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Figure 4. Sensitivity analysis of n in the stochastic model. The size of the perturbation effect, n sitcs , at time f = 20 starting near the disease 
free equilibrium for (A) Between and within-groups infection rates p b and (log scale). (B) Natural mortality rate d. (C) Disease induced mortality rate 
e. (D) Dispersal rate m (log scale). Default parameters are given in Table 2, and one parameter is varied at a time. This is analogous to Fig. 2 for the 
non-spatial case. Initial conditions are such that 20% of sites are randomly chosen to start near the endemic equilibrium (with a minimum of 1 
infective), while the remainder begin at the disease free equilibrium. 
doi:1 0.1 371 /journal.pone.0086563.g004 



state population reduction leads to a persistent perturbation effect, 
as we saw in the non-spatial model (see Fig. S2 in File SI). 
However, this endemic state is very sensitive to the balance 
between site-level establishment and recovery of disease which 
makes it difficult to explore variation in the perturbation effect 
with respect to the value of key parameters. We therefore focus 
attention on the transient perturbation effect when starting close to 
the disease free state in the spatial model. 

Transient perturbation effect in the spatial model 

The behaviour during population reduction in the spatial model 
is shown for the population values S(t), 1(f), and N(t) (see Fig. 3A), 
and the proportion of infected sites Pi(t), dispersal rates and 
effective transmission rate jff eff (see Fig. 3B). The distribution of 
infectives between sites is shown in Fig. 3G before, during, and 
after population reduction. Prior to population reduction, sites can 



be classified as disease-free or infected. During population 
reduction, the typical level of disease within sites decreases, but 
the number of infected sites increases. When population reduction 
ceases, typical prevalence in infected sites returns to previous levels 
which, given that there are now more of them, leads to a rapid 
increase in global prevalence. Some light may be shed on the 
mechanisms behind such changes, as population reduction leads to 
a large increase in dispersal, followed by increasing rates of 
horizontal disease transmission, H (see Fig. 3B). Population 
reduction disrupts the stable demographic structure (shown in 
Fig. 3C) leading to an increase in the dispersal rate and movement 
of infectives to previously disease free sites. This vacuum effect [5,26] 
emerges from the spatial model's density dependent dispersal and 
leads to increased transmission. The effective horizontal transmis- 
sion rate parameter varies with population reduction effort p\. for 
small pi, there is an almost linear increase in P el f (see Fig. 3D), 
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which agrees well with the explicit increase assumed to be 
Pefl = fi + kp in the non-spatial model. However, one key 
difference (as shown in Fig. 3B), is that the increase is not 
immediate, but grows linearly with time — an effect not accounted 
for by our earlier analysis. 

We now explore the sensitivity of this perturbation effect with 
respect to key aspects of demography and disease dynamics. The 
results are broadly consistent with those obtained when starting 
close to the disease free equilibrium in the non-spatial model. Fig. 4 
shows results for parameters analogous to those in Fig. 2, and that 
II decreases with mortality rates d and e and increases with 
dispersal rate m (similar to k in the non-spatial case). The role of 
disease transmission is more complex. The perturbation effect 
decreases with between-group infection rate P b which reduces the 
number of disease free sites, and increases with /?,„, which increases 
disease persistence within sites. Thus for small fij, and sufficient /?„,, 
population reduction is able to spread disease to uninfected sites 
where it is able to persist. We also explore the impact of varying 
the threshold parameter a that determines how sensitive the rate of 
dispersal is to local reductions in the size of the population in the 
destination site (see Appendix S3.1 and Fig. SI in File SI). Results 
show that a perturbation effect occurs for a wide range of values, 
although the largest effects are seen for a around 0.9 (we suspect 
the largest increase would be observed for a near 1 —pi)- 
Perturbation effects were also found for alternative forms of the 
density dependent dispersal function f{Nk) (results not shown). 

Discussion 

In this paper we explored the impact on disease control of 
enhanced transmission resulting from individual or demographic 
responses to population reduction. Using a generic non-spatial and 
deterministic model of demography and disease dynamics we 
explored the potential for such effects to reduce and reverse the 
disease control benefits of population reduction. We found that 
there was a threshold of enhanced transmission above which a 
perturbation effect occurred, whereby the number of infected 
individuals increases during the period when population reduction 
is applied. However, sufficient population reduction (the level 
rising with mortality rates d and e and disease enhancement k, but 
decreasing with infection rate ft) will always reduce numbers of 
infectives in the area it is applied. Disease systems with low levels 
of disease are more sensitive to the impacts of enhanced 
transmission. For systems with endemic disease, the potential for 
the perturbation effect increases with natural and disease induced 
mortality rates (due to reduced levels of endemic disease), with the 
opposite trend where disease is emerging, as higher mortality 
removes cases caused by enhanced transmission. With respect to 
the horizontal transmission rate, the perturbation effect in 
endemic disease is maximised for small to intermediate ft (at the 
point where the disease changes its ability to persist in the absence 
of population reduction). For emerging disease however, higher jS 
causes the disease to reach equilibrium sooner, but reduces the size 
of the perturbation effect, so the earlier the disease is measured, 
the higher the optimal /?, but the weaker the perturbation. 
Enhanced transmission effects can also lead to disease being 
maintained by population reduction in systems where it would 
otherwise die out. 

We also considered a spatially explicit model that represents 
demographic fluctuations and disease transmission within locally 
well mixed populations, and dispersal and disease transmission 



between such groups. In this context we found that enhanced 
transmission emerged implicitly as a demographic response to 
population reduction when dispersal was density dependent. This 
enhancement would be increased if individuals explicitly changed 
their behaviour, e.g. by dispersing more or by increasing agonistic 
interactions and therefore disease contacts between groups (i.e. 
increasing However, the implicit dispersal mechanism alone 
was sufficient to give rise to a perturbation effect. We found that 
the system was susceptible to enhanced transmission in both the 
case of endemic and emerging diseases when infection was 
heterogeneously distributed among groups and when overall levels 
of disease were relatively low. For emerging disease we showed 
that the impact of mortality rates was qualitatively similar to the 
predictions of the non-spatial analysis. In the spatial model, 
dispersal rate played a similar role to the non-spatial enhancement 
parameter k, whereas the role of horizontal disease transmission is 
not directiy comparable between the two cases. In the spatial 
context, higher within-group transmission increased the size of the 
perturbation effect, but even low rates of between-group 
transmission reduced it. Analysis of the effective contact rate in 
the spatial model reveals that enhanced transmission varied in 
time and this could be incorporated in future analysis of the non- 
spatial system. It is worth noting that the linear assumption hp for 
disease enhancement is reasonable, at least early in disease 
intervention period, and for low removal rate p. 

Many authors have noted problems related to disease control 
via population reduction in wildlife [10,21,22,24,32], including 
situations were disease risks are increased rather than reduced 
[13,25]. Individual behavioural [18-20] and demographic [33] 
responses to population reduction are thought to enhance disease 
transmission in wildlife. One system of particular relevance is TB 
in badgers, where the disease is spatially heterogeneous [34], 
transmission between groups is weak, and the host exhibits density 
dependent dispersal [29,30,35]. Thus TB in badgers is a disease 
system exhibiting many of the properties this paper shows are 
likely to lead to the perturbation effect, and we note that the 
RBCT [23] did indeed show that culling was associated with an 
increase in bTB prevalence in badgers [25]. Moreover, the results 
of this paper suggest that a wide range of wildlife disease systems 
are sensitive to such effects, and this is consistent with the marked 
inefficiencies of population reduction as a disease control strategy 
observed to date. However, the effects studied here are likely to be 
even more widespread than current empirical studies suggest as 
they undermine the efficacy of population reduction measures 
even in situations where they do not lead to a complete reversal of 
its effectiveness. 
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